ABSTRACT
The infectious power of several SARS‐COV‐2 variants, commonly described by R0, are rationalized under a molecular interaction perspective using molecular dynamics simulations, quantum mechanical descriptors of the nature of bonding interactions, and calculated spectral features of the specific mutations in the spike protein and the interactions of the replaced residues with the ACE2 receptor. More information can be found in the Full Paper by C. Cappelli, A. Restrepo et al.
ABSTRACT
Specific S477N, N501Y, K417N, K417T, E484K mutations in the receptor binding domain (RBD) of the spike protein in the wild type SARS-COV-2 virus have resulted, among others, in the following variants: B.1.160 (20A or EU2, first reported in continental Europe), B1.1.7 (α or 20I501Y.V1, first reported in the United Kingdom), B.1.351 (ß or 20H/501Y.V2, first reported in South Africa), B.1.1.28.1 (γ or P.1 or 20J/501Y.V3, first reported in Brazil), and B.1.1.28.2 (ζ, or P.2 or 20B/S484K, also first reported in Brazil). From the analysis of a set of bonding descriptors firmly rooted in the formalism of quantum mechanics, including Natural Bond Orbitals (NBO), Quantum Theory of Atoms In Molecules (QTAIM) and highly correlated energies within the Domain Based Local Pair Natural Orbital Coupled Cluster Method (DLPNO-CCSD(T)), and from a set of computed electronic spectral patterns with environmental effects, we show that the new variants improve their ability to recognize available sites to either hydrogen bond or to form salt bridges with residues in the ACE2 receptor of the host cells. This results in significantly improved initial virusâ â â cell molecular recognition and attachment at the microscopic level, which trigger the infectious cycle.
Subject(s)
COVID-19 , SARS-CoV-2 , Humans , Mutation , Protein Binding/genetics , SARS-CoV-2/genetics , Spike Glycoprotein, Coronavirus/chemistryABSTRACT
The magnified infectious power of the SARS-CoV-2 virus compared to its precursor SARS-CoV is intimately linked to an enhanced ability in the mutated virus to find available hydrogen-bond sites in the host cells. This characteristic is acquired during virus evolution because of the selective pressure exerted at the molecular level. We pinpoint the specific residue (in the virus) to residue (in the cell) contacts during the initial recognition and binding and show that the virusâ â â cell interaction is mainly due to an extensive network of hydrogen bonds and to a large surface of noncovalent interactions. In addition to the formal quantum characterization of bonding interactions, computation of absorption spectra for the specific virusâ â â cell interacting residues yields significant shifts of Δλmax =47 and 66â nm in the wavelength for maximum absorption in the complex with respect to the isolated host and virus, respectively.
Subject(s)
Angiotensin-Converting Enzyme 2/metabolism , COVID-19/prevention & control , Receptors, Virus/metabolism , SARS-CoV-2/metabolism , Spike Glycoprotein, Coronavirus/metabolism , Amino Acid Sequence , Angiotensin-Converting Enzyme 2/chemistry , Angiotensin-Converting Enzyme 2/genetics , COVID-19/epidemiology , COVID-19/virology , Humans , Molecular Dynamics Simulation , Pandemics , Protein Binding , Protein Domains , Receptors, Virus/chemistry , Receptors, Virus/genetics , SARS-CoV-2/genetics , SARS-CoV-2/physiology , Spike Glycoprotein, Coronavirus/chemistry , Spike Glycoprotein, Coronavirus/geneticsABSTRACT
We present a detailed computational study of the UV/Vis spectra of four relevant flavonoids in aqueous solution, namely luteolin, kaempferol, quercetin, and myricetin. The absorption spectra are simulated by exploiting a fully polarizable quantum mechanical (QM)/molecular mechanics (MM) model, based on the fluctuating charge (FQ) force field. Such a model is coupled with configurational sampling obtained by performing classical molecular dynamics (MD) simulations. The calculated QM/FQ spectra are compared with the experiments. We show that an accurate reproduction of the UV/Vis spectra of the selected flavonoids can be obtained by appropriately taking into account the role of configurational sampling, polarization, and hydrogen bonding interactions.